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Abstract 

We discuss the stability of multiquark systems within the recent 
model of Glozman et al. where the chromomagnetic hyperfme in- 
teraction is replaced by pseudoscalar-meson exchange contribu- 
tions. We find that such an interaction binds a heavy tetraquark 
systems QQqq (Q = c, b and q = u, d) by 0.2 — 0.4 GeV. This is 
at variance with results of previous models where ccqq is unsta- 
ble. 
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The existence of tetraquark hadrons — two quarks and two antiquarks — has been raised 
about twenty years ago by Jaffe [1] and has been studied within a variety of models. The 
MIT bag study indicated the presence of a dense spectrum of tetraquark states in the 
light sector [1]. Later on, tetraquark systems have been examined in potential models [2-4] 
and flux tube models [5]. In particular the question of stability has been raised, that is 
whether the tetraquark ground state lies below or above the lowest (qq) + (qq) threshold. 
Weinstein and Isgur showed [2] that there are only a few weakly-bound states of resonant 
meson-meson structure in the light (u, d, s) sector. On the other hand, no bound state 
was found by Carlson and Pandharipande [5] in their flux-tube model with quarks of equal 
masses. 

One important result in the MIT bag or potential models is that the chromomagnetic 
interaction plays a crucial role [2,6] in lowering the ground state energy of a light system. 
Otherwise, in a system of two heavy quarks and two light antiquarks QQqq (Q = c or 
b, q = u, d or s) stability can be achieved without spin-spin interaction, provided the 
mass ratio m(Q)/m(q) is larger [3] than about 15, which means that Q must be at least a 
6-quark. 

In the light sector, there are several candidates for non-qq states, but the experimental 
situation is not yet conclusive. For a review, see for example Ref. [7] and the last issue 
of Review of Particle Properties [8]. In the heavy sector, experiments are being planned 
at Fermilab and CERN, to search for new hadrons and in particular for doubly charmed 
tetraquarks [9-11]. 

Recently, the baryon spectrum has been analysed by Glozman et al. [12,13] within a 
chiral potential model which includes meson-exchange forces between quarks and entirely 
neglects the chromomagnetic interaction. In view of its intriguing success in the description 
of the baryon spectrum, it seems to us natural to apply this model to multiquark hadrons 
with more than three quarks, and in particular to tetraquarks. 

A general Hamiltonian containing both chromomagnetic interaction and meson-ex- 
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change contributions has the form 



- Z ^ • \ c ^ • ^ ^ ) - Z • ^ F *i • ^ ^ ( r ^- ) > 

where r»i is the constituent mass of the quark located at f^; = |Fj — Fj| denotes the 
interquark distance; al, A°, Af are the spin, colour and flavour operators, respectively. 
Spin-orbit and tensor components may supplement the above spin-spin forces for studying 
orbital excitations. The potential in H has three parts containing the confining, the 
chromomagnetic and the meson-exchange contributions, respectively. 

The confining term V^ on f usually consists of a Coulomb plus a linear term, 

V con{ = -- + br + c. (2) 

r 

In the following, we shall use either the very weak linear potential of Glozman et al. [13] 
corresponding to 

(d) a = c = 0, and b = 0.01839 GeV 2 , (3) 
or a more conventional choice 

(C 2 ) a = 0.5203, b = 0.1857 GeV 2 , c = -0.9135 GeV, (4) 

which has already been used in the study of tetraquarks by Silvestre-Brac and Semay [4] . 

The third term in H is often understood as the chromomagnetic analogue of the Breit- 
Fermi term of atomic physics. For mesons, A^-A^j = —16/3, and a positive V g , as in the 
one-gluon-exchange model, shifts each vector meson above its pseudoscalar partner, for 
instance D* > D in the charm sector. For baryons, where A^-A^ = —8/3 for each quark 
pair, such a positive V g pushes the spin 3/2 ground states up, and the spin 1/2 down, for 
instance A > N. In Ref . [4] , the following radial shape has been used 

mimjd z r 
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with the same value of a as in Eq. (4) and d = 0.454 GeV 1 . 



The last term of H corresponds to meson exchange, and an explicit sum over F is 
understood. If the system contains light quarks only (as in Refs. [12,13]), the sum over 
F runs from to 8 (1 — 3^7r, 4 — 7 K, 8 ^ t] and — > rj'). If a heavy flavour is 
incorporated, a phenomenological extension from SU(3)f to SU(4)p would further extend 
the sum to F = 9 — 12 corresponding to a D-exchange, F = 13 — 14 to a D s -exchange and 
F = 15 to an ?y c -exchange. If Vf = 0, one recovers a standard constituent quark model. 
The radial form of Vf ^ is derived from the usual pion-exchange potential which contains 
a long-range part and a short-range one 



(6) 
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where n is the pion mass. A coupling constant g 2 /At: = 0.67 at the quark level corresponds 
to the usual strength g^NN/An — 14 for the Yukawa tail of the nucleon-nucleon (NN) 
potential. 

When constructing NN forces from meson exchanges, one disregards the short-range 
term in Eq. (6), for it is hidden by the hard core, and anyhow the potential in that region 
is parameterized empirically. Similarly, when Tornqvist [14], Manohar and Wise [15] or 
Ericson and Karl [16] considered pion exchange in multiquark states, they had in mind the 
Yukawa term exp(— \ir)jr acting between two well-separated quark clusters. For similar 
reasons, Weber et al. [17], in their model with hyperfine plus pion-exchange interaction 
ignored the delta-term too. Therefore it is somewhat of a surprise to see the delta-term 
of Eq. (6) taken seriously, and with an ad-hoc regularisation playing a crucial role in the 
quark dynamics [12,13]. This regularised form is [13] 

V, = 0(r - ro)/J 2 ex P (-H _ 4^ exp( _ £2(r _ ^ (?) 



with the Yukawa-type part cut off for r < tq, where tq = 2.18 GeV 1 , e = 0.573 GeV, and 
H = 0.139 GeV for n, 0.547 GeV for 77 and 0.958 GeV for rf. 

Incorporating in H both light quarks and charm (one can add bottom similarly) and 
working out the flavour matrix elements of the mesons-exchange interaction between two 
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quarks (or antiquarks), one obtains: 



Kr + -V UU + -V 

3 77 6 
2V K - ^V^ + 2V l 



[2]f I 



1 



+ ; [2] F * 



1/2 



4 3 

_l _ycc 

3 7, 2 77 c 



[2]f I 







(8) 



-2V K - 2 -V^ s - 2VJ 




[11] F /= 



[H]f /= 







1/2 



for F = 1, . . . 15, and 



( v ij) = o°i-°j v v' 



(9) 



for F = 0. This is an extension of Eq. (3.3) of Ref. [12] from SU(3) to SU(4). Here I is 
the isospin and in each case it is specified whether the pair is in a symmetric [2]p state or 
in an antisymmetric [11]f state. Actually little uu or dd mixing is expected in r\ c so that 
the contribution of can be safely neglected. Moreover when the meson mass \x reaches 
values of a few GeV as for D or r\ c the two terms in Eq. (6) basically cancel each other 
so the expressions (7) reduce practically to their SU(3) form [12]. This is in agreement 
with Ref. [18] where it has been explicitly shown that the dominant contribution to the S c 
and E* masses is due to meson exchange between light quarks and the contribution of the 
matrix elements with the D (D s ) and D* (D*) quantum numbers (which are evaluated 
phenomenologically, fitting the mass difference E c — A c ) play a minor role. In the following 
numerical calculations we will neglect the exchange of heavy mesons. 

The remaining parameters are the quark masses. They are indicated in Table 1, in 
conjunction with the two choices. The light quark masses m = m u = ma are from the 
corresponding literature [13,4]. The heavy quark mass tuq = m c or mj is adjusted to 
reproduce the experimental average mass M = (M + 3M*)/4 of M = D or B mesons. 
The meson mass is obtained from a trial wave function of type <p exp(— ar 2 /2), with a 
as a variational parameter. It has been checked that the error never exceeds a few MeV 
with respect to the exact value. The variational approximation is retained for consistency 
with the treatment of 3- and 4-body systems discussed below. 
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Table 1. Quark masses and average heavy meson masses M = 
(M + 3M*)/4 (M = D, B) used for the potentials (d) and 
(C 2 ). The units are GeV. 



Model 


m 


m c 


nib 


D 


B 


(Ci) 


0.340 a ) 


1.350 


4.660 


2.001 


5.302 


(C 2 ) 


0.337 b ) 


1.870 6 ) 


5.259 b) 


2.006 


5.350 



■) Ref. [13] ^ Ref. [4] 



We now briefly discuss the baryons. In the model of Glozman et al. the explicit form 
of the Hamiltonian integrated in the spin-flavour space is : 

2 f 15V 7r -V r) -2(g /g) 2 V r] , for N 
H = H + -^{ (10) 
48 ^ 2 I 3V 7V + V ri + 2(g /g) 2 V v , for A 

with 

i i<j 

where g 2 /4n = 0.67, (g /g) 2 = 1.8 and = V^ u of Eq. (8). 

We have performed variational estimates with a wave function <p oc exp(— a(p 2 +A 2 )/2), 
where p = r 2 — F3 , A = (2?i — r 2 — F3) / v^J, and reproduced the results of the more elaborate 
Faddeev calculations of Ref. [13]. When the meson-exchange terms are switched off, the iV 
and A ground states are degenerate at 1.63 GeV. When the coupling is introduced, the wave 
function is modified. For the nucleon, the spin-independent part Hq of the Hamiltonian 
gives a contribution of 2.11 GeV, and receives a large —1.14 GeV correction from meson 
exchange. For the A ground state, the contribution of Hq and meson exchange parts are 
1.91 GeV and —0.63 GeV, respectively. Thus one ends up with a reasonable value for the 
A- N splitting, close to 0.3 GeV. 

We have also calculated the ground state baryons of content cqq with a trial wave 
function <p oc exp(-(ap 2 + /3A 2 )/2) and found A c = 2.32 GeV and E c = E* = 2.48 GeV, 
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close to the experimental values and consistent with the findings of Ref. [18], although the 
Hamiltonian, its treatment, and the input parameters are somewhat different there. 



Due to arguments at the beginning of this letter, here we discuss tetraquarks contain- 
ing heavy flavours, i.e. QQqq, studying the most favourable configuration 33, S = 1, I = 0. 
This means QQ is in a 3 colour state and qq in a 3 colour state. The mixing with 66 is 
neglected because one expects this plays a negligible role in deeply bound heavy systems 
[3]. Then the Pauli principle requires S12 = 1 for QQ, and S34 = 0, I34 = for qq, if the 
relative angular momenta are zero for both subsystems. This gives a state of total spin 
S = 1 and isospin 1 = 0. 

The tetraquark Hamiltonian integrated in the colour-spin-flavour space, and incor- 
porating the approximations discussed in relation to Eq. (8), reduces to 

H = 2(m + m Q ) + ^ + ^ + ^±^-f z + J2Vi j , (12) 
rriQ m zmrriQ ^ 



where 







( _ a 


Vi2 = 










V r 12 






' a 


Vy = 


l( 








( a 


V34 = 


3l 








V r 34 



+ 6ri2 + c 
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i = 1 or 2, j = 3 or 4, (13) 



The momenta p x , etc., are conjugate to the relative distances x = r*i — F2, y = F3 — F4, and 
z = (?! + r 2 - r 3 - r 4 )/v / 2- The wave function is parameterized as 

^ oc exp[-(ax 2 + (3y 2 + 7 ^ 2 )/2], (14) 

and the minimization with respect to a, (3 and 7 leads to the results displayed in Table 2 
for Q = c and b. Columns 2 and 3 corresponds to results derived from Eqs. (12) and (13). 
This shows that both the ccqq and bbqq systems are bound whatever is the potential, (Ci) 
or (C2), provided meson exchange is incorporated. This is in contradistinction to previous 
studies based on conventional models where the flavour-independent confining potential is 
supplemented by one gluon exchange. This remark is illustrated by column 4 which shows 
that ccqq is unbound in such a case. 
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Table 2. Energy (GeV) AE = QQqq - 2(Qq) for heavy tetra- 
quarks with Q = c or 6, q = u, d in three different models with 
parameters defined in the text (OME = one meson exchange). 



System 


(Ci)+ OME 


(C 2 )+ OME 


Ref. [4] 


ccqq 


-0.185 


-0.332 


0.019 


bbqq 


-0.226 


-0.497 


-0.135 



Considering this result one could rise the question if in the model of Glozman et al. 
a proliferation of multiquarks systems appears. We have therefore tried to investigate 
QQqqqq and g 6 systems as well. We have proceeded analogously to the previous case using 
a Gaussian wave function. In the case of QQqqqq, group-theory analysis [19] shows that 
the most favourable configuration is the one where the light quark subsystem has S = 1, 
1 = 0, corresponding to a global spin-flavour-averaged interaction {V} = lOV^ — 2/3 — 
4/3(go/g) 2 V rj ' . Our numerical calculation shows that this potential is largely insufficient to 
bind the system. Also the q 6 system, the most favourable configuration of which is S = 1, 
7 = 0, leading to (V) = HV n — 5/3V v — 10/3(go/g) 2 V ri ' , is not sufficiently bound for being 
under the two baryons threshold. 

It is also interesting to notice that the binding energy of ccqq and bbqq are nearly twice 
larger for the potential (C2) as compared to those of {C\) and also much more different from 
each other. The reason is that (C2) contains a Coulomb part which binds more, heavier 
is the system, leading thus to a larger separation among levels as well. The potential (C 2 ) 
has been fitted to reproduce the J/^> and the T meson masses. It also gives overall good 
results both for other mesons and baryons. By construction [12,13], the potential (Ci) was 
designed and fitted to light baryons only. It is desirable to construct a chiral potential 
model with a wider range of validity, covering the light and heavy sector as well, and to 
apply it to the study of mesons and multiquarks systems. Our results indicate that the 
chiral model of Glozman et al. leads to qualitatively different results for tetraquarks, with 
respect to commonly used quark models. We hope that future experimental investigations 
might distinguish among various approaches to quark dynamics. 
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